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Abstract 

We investigate the prospect of stop squark search at Tevatron RUN-II in 
mSUGRA motivated as well as unconstrained supersymmetric models, when 
the lighter stop squark turns out to be the next lightest supersymmetric par¬ 
ticle (NLSP). I 11 this case the decay into a 4-body final state consisting of a b 
quark, the lightest neutralino and two light fermions may compete with the much 
publicized loop induced two body decay into a charm quark and the lightest neu¬ 
tralino. We systematically study the parameter space in mSUGRA where the 
lighter stop squark turns out to be the NLSP and calculate the branching ratios 
of the competing channels in both models. Our results show that the four body 
decay may indeed be the main discovery channel particularly in the low tan (3 
scenarios. We discuss the detectability of stop squark pairs in the 4-body decay 
channel leading to one lepton with 2 or more jets accompanied by a large amount 
of missing energy. We also studied the corresponding background processes and 
the kinematic cuts required to suppress them using parton level Monte Carlo 
simulations. We have commented upon with illustrative examples, the required 
revision of the existing mass limits of the stop NLSP assumed to decay solely 
into the loop induced 2-body channel in the presence of the competing 4-body 
decay. 
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I. Introduction 


The Minimal Supersymmetric Standard Model (MSSM) |l|] is a well motivated extension 
of the Standard Model(SM). As of now, neither there is any evidence of it nor has it been 
ruled out by the electroweak precision measurements at LEP ||. Unfortunately, we are not 
equipped with any information about the range of superparticle masses from the theoretical 
point of view. On the other hand, from the experimental side, there are some lower bounds 
from the non observation of superparticles in colliders, such as LEP || and Tevatron RUN- 

i0- 

The up type squark of the third generation - the stop squark, the superpartner of the top 
quark - is something special. It is because, of the large top Yukawa coupling which controls 
the evolution of the soft-supersymmetry breaking masses of the left and right handed stop 
squarks, £/,,£#, via the Renormalisation Group equations. This tends to reduce these masses , 
compared to the other squark masses 0. Moreover, because of the large top quark mass, the 
two weak states £r may mix very strongly leading to a relatively large splitting between the 
two physical mass eigenstates fi ,^2 @ (in our notation m > m t -J. Interestingly, the mass 
of the lighter states t\ may be even below the top mass. In fact, it is quite conceivable that 
it happens to be the next lightest supersymmetric particle (NLSP), the lightest neutralino 
Xi being the lightest supersymmetric particle (LSP) by assumption. 

In the coming years the second phase of Tevatron experiments, the RUN-II, will start 
operating with an integrated luminosity of at least 2 fb -1 per experiment at 2 TeV center 
of mass energy, which is about ten times larger than the acquired luminosity in RUN-I with 
center of mass energy 1.8 TeV. With a further luminosity upgrade it is expected to collect a 
luminosity of something like 15-20 fb _1 after a few years of operation. The phenomenology 
of the stop squark could be of special interest, since it might be the only strongly interacting 
sparticle within the kinematic range of RUN-II experiments. It is, therefore, very important 
to fix up the strategies for isolating the stop signal for all conceivable decay modes. 

Currently, the search for stop at LEP [||] and Tevatron RUN-I 0 experiments have 
yielded negative results. The most stringent bound comes from Tevatron experiments which 
puts a lower limit on lighter stop mass m- tl >119 GeV for o =40 GeV and this limit 
becomes slightly weaker for higher value of rri^o , e.g, rnq > 102 GeV for m^o =50 GeV [0. 
In deriving these limits, it was assumed that the loop induced, flavor changing decay into a 
charm quark and the LSP 

U -> CXI (1) 

occurs with 100% branching ratio(BR) - until recently assumed to be a valid assumption 
if the t\ state happens to be the NLSP. In the R parity conserving model, LSP cannot 
decay further and escapes the detector resulting a large imbalance of transverse energy. The 
stop signal is tagged through identifying jets and missing energy. Since the production of 
stop pairs is dominantly via QCD and depends on its mass only, the above limits from the 
Tevatron are fairly model independent, except for the dependence on rn^o , which influences 
the efficiency of the kinematical cuts. 
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In the minimal supersymmetric extension of the standard model with arbitrary soft break¬ 
ing masses of the superparticles, the stop may be the NLSP. However, the stop NLSP may 
also be realized in more constrained models like the minimal version of the supergravity 
model (mSUGRA). The latter is the most economical model in the context of supersymme¬ 
try (SUSY) searches in colliders. 

In the mSUGRA model, the supersymmetry breaking takes place in hidden sector which 
is communicated to the visible sector by gravitational interactions. The mass spectrum of 
sparticles at a lower energy scale can be obtained by the Renormalisation Group equations 
from the inputs at some higher scale, which is usually assumed to be the Grand Unification 
theory(GUT) scale (Mg). In mSUGRA models, these input at Mg are the common scalar 
mass (mo), the common gaugino mass (mi/ 2 ), the tri-linear scalar coupling term (Ao), tan/? 
(the ratio of two vacuum expectation values of the two higgs doublets which generate the 
masses of the fermions and gauge bosons through electroweak symmetry breaking) and the 
sign of /i, (the higgsino mass parameter) 0. The magnitude of fi is fixed by the radiative 
electroweak symmetry breaking condition (REWSB) ||. Surprisingly due to the interplay 
of these parameters, which we will discuss later, there is a substantial region of mSUGRA 
parameter space where the t\ happens to be the NLSP and sometimes the only squark within 
the striking range of the Tevatron. 

The authors of 


10, 11] have emphasized that a competiting channel of t\ decay may be 


there even if the stop is the NLSP. The 4-body decay of t\ into a b quark, the LSP and two 
approximately massless fermions ||T0], |TTf, 




( 2 ) 


via heavier SUSY particles, may also have significant decay rates, or even dominate over 
the loop induced decay mode, eq.[l], for certain regions of the MSSM or mSUGRA parameter 
space. As a consequence, this 4-body decay mode along with stop pair production in colliders 
may yield final states containing a lepton plus jets accompanied by huge amount of missing 
energy Q This event topology appears to be indeed different in comparison to the final state 
resulting from the loop induced stop decay eq. |T]. The new observation on the 4-body decay, 
therefore, necessitates revision of the strategy to detect the signal of t\ when it is the NLSP. 

In the framework of most of the SUSY models, stop squarks has many interesting decay 
modes depending on its mass. If it is sufficiently heavy, then the main decay modes occurs 
through top quarks and neutralinos, 

U -»• tx°j (j = 1 - 4) (3) 

if it is kinematically allowed. There is also another competitive decay channel through 
charged current interactions into bottom quarks and a lighter chargino(xf), 

h -> bxt (4) 

4 It is to be noted that similar final states occur from the production of other SUSY particle ( e.g., squarks, 
gluinos) and their subsequent cascade decays in colliders Q. These signals are used in conventional SUSY 
searches. 
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If these modes are not kinematically accessible, but still heavy enough, then the 3-body 
decay modes to bottom quarks, a W boson or a charged higgs scalar H ± and neutralinos, t\ —> 
bW + x°j or ti —> bH + x > j (j=l-4), can be accessible [|I2[] . These decay channels, particularly, the 
W boson final states may be dominant in absence of the 2-body decay channels, eq.(|3|{4|). 
Beside these 3-body decay modes, in the light slepton scenario which may arise in some 
mSUGRA models, lighter stop can decay with a final states containing sleptons H, [TO, O, [14 


t\ —> Mu ; ti —> Mu 


( 5 ) 


The search prospects of ti state at Tevatron experiments has been investigated by many 
authors [|], [L3|, |l^, |l^. [17]. [R|. These studies were carried out by considering the 2-body 
decay mode of i\ into a b quark and a lighter chargino, eq.[|, yielding a single lepton or a 
dilcpton pair plus large amount of missing energy with some hadronic activities following the 
cascade decays of xt info a LSP and massless fermions, xt XifP■ Though this final state 
has the same particle content as in eq.2, the kinematical characteristics are quite different 
in the two cases. Hence a slightly different set of kinematical cuts is needed for isolating the 
channel, in eq.2, from the background. 

Similar studies were also carried out for the same signal with different kinematics in 
the light slepton scenario where ti decays via 3-body into a b quark and a charged slepton 
(lepton) and neutrino (sneutrino), eq.|5|. The SUSY searches take a dramatic turn in the high 
tan (3 regime where lighter staus(fi), the SUSY partner of tau lepton turns out to be very 
light, even may be the NLSP. It yields a huge number of tau leptons in the SUSY particle 
production and their cascade decays in colliders [|T9|. In Ref. [18|] discovery potential of t\ 
states has been studied for Tevatron RUN-II in high tan /3 regime. 

However, the discovery potential of t\ of the stop NLSP decaying dominantly into the 
4-body final state, especially in the context of the mSUGRA model, has not yet been studied 
systematically for the upgraded Tevatron. The purpose of this work is to investigate this 
possibility further. In our study, we concentrated on the channel with a single lepton plus 2 
or more jets accompanied by large amount of missing energy. This final state appears from 
stop pair production following the semileptonic decay of one t\ and hadronic decay of the 
other. We also discuss the possible sources of SM backgrounds and optimize the set of cuts 
to minimize the contamination. 

Our investigation is based on models which conserves R parity. In R parity breaking 
models, few more channels of t\ decay open up which depend on R parity violating couplings. 
The collider phenomenology of t\ in presence of R parity breaking SUSY model will be 
discussed elsewhere [2H| . 

We have organized our paper as follows. In Sec. II, we have isolated the region of the 
parameter space where lighter stop happens to be the NLSP in the mSUGRA model. In Sec. 
III, the key points of stop decay patterns have been mentioned. The parameter spaces, both 
in mSUGRA and MSSM, where the branching ratio of the 4-body decay is significant, are 
identified. In Sec. IV, we investigate the signal of t\ in the 4-body decay channel by triggering 
single lepton plus jet with missing energy channel. We also estimate the corresponding SM 
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backgrounds after minimizing them by imposing kinematical cuts. We conclude in Sec.V 
with a few remarks. We have also commented on the possibility of revision in the current 
mass limits in the presence of the competeting 4-body decay mode along with illustrative 
examples. 


II. Stop NLSP in the mSUGRA model 


In this section we focus our attention on the parameter space in the mSUGRA model where 
the stop is the NLSP. As discussed in the introduction the model is chracterised by five free 
parameters: m 0 , mi/ 2 , A 0 , tan/3 and sign(/r). 

For a given mo, mi/ 2 , tan/3 and sign(/r), the stop turns out to be the NLSP for a range 
of A 0 (| Aomin | < | A)[ < | Ao max |). The sign of A 0 is required to be the opposite to that of fi 
in order to produce larger mixing in the stop mass matrix. For |A 0 | < |A 0 mm| the mixing is 
small and the lighter stop becomes heavier than the chargino. For |Ao| > |Amax|, the lighter 
stop becomes lighter than the lightest neutralino which is also the lightest supersymmetric 
particle. An additional constraint on |A 0ma x| exist in principle due to the requirement that 
the vacuum must not break charge and color symmetry (the CCB condition) (21J. However, 


quite often this constraint becomes redundant in the mSUGRA model since the requirement 
that the stop be heavier than the LSP produces a stronger bound at least in the low mo, 
m 1/2 scenario. Moreover the upper bound is in general unimportant for this paper, since the 
BR of the 4-body decay is significant only when mi is close to m-± so that the chargino is 
of relatively small virtuality. This happens for A 0 close to A 0min . 

For illustration, in Fig.l we present the stop mass as a function of A 0 for different tan/3, 
m 0 =200, mi/ 2 = 145 , and /i > 0 (all masses and mass parameters in this paper are in 
GeV). Here 777 -±=106 for tan/3=5 which varies modestly with tan/3. Only those values of 
A are considered for which the t\ is the NLSP. As seen from the figure A 0rnax nearly equal 
to 565 and \A 0m i n \ nearly equal to 510 for tan/3=3. For a given A 0 , the mass of the lighter 
stop is larger for larger tan/3 since mixing is reduced by the in the stop mixing matrix. 

Similar information for t? 7 o= 140 , m 1 / 2 =180, fi > 0 is given in Fig.2. In Fig.2 we find 
that higher masses of the stop NLSP are allowed since here is somewhat larger. Here 

??7-±=136 for tan/3=5 which varies modestly with tan/3. 

If ?77 0 >> mi/2i then mq tends to increase, unless A 0 is properly tuned so that a large 
mixing in the stop mass matrix still leads to a low mass stop NLSP. Beyond some large value 
of ?77 0 extreme fine tuning is needed which is rather unaesthetic. 

On the other hand for mi/ 2 » m 0 , quite often a large A 0 is needed to produce a stop 
NLSP, which is in conflict with the invariance of the vacuum under charge-color symmetry 
21], Even if the CCB condition is disregarded, we find that in a very narrow range of 


Ao we have the stop NLSP. Otherwise a slepton NLSP scenario is obtained over most of 
the parameter space. Thus in this paper we shall focus our attention on scenarios with 
777 0 ~ ?77i/ 2 . 

It is clear from Fig.l and Fig.2 that if the BR of the loop induced decay is indeed ~ 
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100% as is normally assumed, then the limits from LEP and Tevatron are already sufficiently 
restrictive and can exclude a large region of the mSUGRA parameter space ( in particular 
certain ranges of Aq for given values of other mSUGRA parameters will be excluded ). It is 
also reasonable to hope that RUN-II at a much higher luminosity will extend this probe to 
a much larger region of the parameter space. However, as emphasized in the introduction, 
the 4-body decay, eq.2, can indeed compete with the loop induced 2-body decay and reduce 
its BR significantly [?, |TT[. This will be illustrated with several examples in the context of 
the mSUGRA model in the next section. The existing limits on rrij A , therefore, may require 
a revision. We shall try to illustrate estimate the required revision in section III. 
More importantly, when new and stronger limits come from RUN-II, the interplay of the two 
competing channels should be kept in mind. Such limits will depend not only on rn^ and 
m^o, but also on the relative BRs of the two competing channels for stop decay and, hence, 
on other SUSY parameters. 


III. Stop decay branching ratios 


As mentioned in Sec.I, when the t\ is the NLSP, two decay channels are allowed which are 
competitive with each other - the loop induced flavor changing decay mode, eq.[lj 0 and 4- 
body final states with nearly massless fermions, the bottom quark and the LSP, eq. 10.0. 

The 4-body decay mode of lighter stop, eq.^j occurs through many diagrams mediated by: 
(a) W and/or H ± with virtual t, b and X 121 (b) virtual xf 2 with £ or zy. The dependence of 
4-body decay rates on supersymmetric parameters has been discussed in great detail in the 
paper of Ref. |[T0| , |TT] . Nevertheless, few comments are in order: 


• Among all the diagrams, only few have significant contributions. As for example, 
the diagrams mediated by H ± and W are heavily suppressed because, m H ± are larger 
than M w and Yukawa coupling between H ± and fermions are suppressed by respective 
fermion masses. In addition to these, the diagrams mediated by squarks have very little 
contribution as those are expected to be larger(>250) |^2| than t\ state. The top quark 
mediated diagrams do not give large contributions unless ~m t + m^o. 


The diagrams which expected to give a significant contribution to the 4-body decay rate 
are mediated by charginos xf (i= 1,2) and sleptons (£, uf). When the virtuality between 
t\ and xf , £ or is not very large then these diagrams do give a large contributions. As 
for example, the diagrams involving tau sleptons which could be rather light because 
of mixing (especially, for high tan/3 |18|) may contribute appreciably. 


In the MSSM the soft breaking parameters can be chosen arbitrarily. One can, therefore, 
suppress the loop decay to the extent one requires. Now if the stop is chosen to be the 
NLSP, the 4-body decay will automatically have a large branching ratio. For the detectable 
signal, which we consider, one also needs a significant leptonic BR of the stop (see section 
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IV ). This can be achieved by choosing relatively light sleptons. In the following we shall 
illustrate this with different inputs. 

What, however, is interesting is that even in a constrained model like mSUGRA, one 
can have a significant BR for the 4-body decay in some appropriate region of the parameter 
space. Moreover, a significant leptonic BR can be accommodated in a small but interesting 
region. 

We first consider the unconstrained MSSM model. As is well-known the loop decay width 
is controlled by the parameter e which denotes the amount of £l,r _ Cl mixing || and enters 
in the decay width formula as, 

T(ti ex?) = j\e\ 2 fm h ^1 - -jf'j (6) 


The detailed expressions for e and the function / can be found in Ref. j[|. [TTj and Ref. [JTO], [13 
respectively. If the mixing angle 9 t in the stop sector is chosen appropriately so that it 
yields very small value of the e parameter, the loop decay is reduced drastically [[10]. [TT . 
This effect is shown in Fig.3. The choice of the SUSY parameters for this figure are as 
follows assuming gaugino mass unification: the SU(2) gaugino mass parameter M2=180, the 
higgsino mass parameter // = +400, t an/3=4. This three parameters completely defines the 
Chargino and Neutralino sector. We have taken CP-odd neutral higgs mass M ^=150 which 
requires to calculate the neutral higgs to obtain the value of the e parameter. The others 
MSSM parameters requires for calculating the stop branching ratios are the following: the 
common scalar squark mass m ? -=500, common slepton mass m^=175, the trilincar coupling 
in the bottom sector A;,=300, the trilinear coupling in the tan sector A T =200. With these 
parameters, cos 6^=0.11 which leads to a very small e. Because of this the loop decay BR 
sharply falls as soon the 4-body decay is kinematically allowed. The total BR for the four 
body decay exceeds the loop decay BR for rrqq ~ 95 and approaches 100% for slightly higher 
values of m^. For our parameter set m~± rs 162. Thus the 4-body BR can be rather large 
even for charginos of relatively high virtuality. We have checked that even if this smaller 
value of the angle Q- t is not used so that e is somewhat larger, the enhancement of the 4-body 
decay BR, though less dramatic, still exist over a significant region of the parameter space. 
Moreover the leptonic BR of the stop is greater than 10%. This, as we shall see in the next 
section, as adequate for the detection of the 4-body decay of the stop in the channel studied 
by us. 

We are now in a position to rexamine the possible impact of the large BR of the 4-body 
decay channel on the current and future limits on . For the choice of parameters in Fig.3, 
777^0=86.41. For such a large value of the LSP mass, no limit on is currently available. 
However, the expected mass limits from at RUN-II needs careful handling. As we see from 
Fig.3 the total BR of the 4-body decay approaches 100% for rn^ > 100. Obviously, the 
searches for t\ NLSP in this mass range should be based on the 4-body decay channel. 

We now turn our attention to the 4-body BR in the mSUGRA model. We present the 
4-body BR as a function of (Fig. 4 and 5). The choices of SUSY parameters are as 
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Figure 3: The branching ratio BR(ty —> bxiff ') and BR(U —>■ cy?) as a function of the stop squark 
mass (mf ) in MSSM scenario, where the common scalar squark mass m,j=500, common slepton 
mass mf= 175, M 2 =180, p = 400, ^=300, R r =200, M^=150. tan/?=4, and cosg t = 0.11. 


in Fig.l and Fig.2 respectively. In Fig.4, m-± ~ 103, although its precise value depends 
on tan/?. We find that the BR is appreciable when low values of tan/? are considered and 
the virtuality of the chargino is small. This is because the loop decay rapidly increases with 
tan/?. Thus total 4-body BR > 10% is allowed only for tan/? ~ 7. Thus the phenomenology of 
the stop NLSP and its 4-body decay is interesting for a rather limited region of the H 0 -tan/? 
parameter space. 

For the set of parameters in Fig.4, m^o ^55. In this scenario the current mass limit 
from the Tevatron based on the loop decay is rn^ -100. However, in the tan/?=3 scenario 
the total BR of 4-body decay is quite sizable for -85. Hence masses significantly lower 
than 100 cannot be excluded apriority. It follows from Fig.4 that for tan/? - 5 the 4-body 
decay may indeed play a significant role in the discovery of ty. In fact for the entire range 3 
- tan/? - 5 both the decay modes eq. (B® of the stop NLSP should be taken into account. 
Negative results from stop search on the other hand, already excludes certain ranges of the 
A 0 parameter. For example, the low BR of the 4-body decay for tan/? >5 indicate that the 
current limits from Tevatron RUN-I based on the loop decay of ty are valid for these tan/?. 
Referring back to Fig.l, we find that, e.g, A 0 -525 is excluded for tan/?=7 and mo, mi/ 2 , 
sign(/i) as given above. This limit is certainty stronger than the one obtained from the CCB 
condition or from the requirement that stop be heavier than the LSP. Moreover, the limit 
from the CCB condition does not hold if the SU(2)<g)U(l) breaking minimum happens to be 
a false vacuum with a life time larger than the age of the Universe [23]. In contrast the limits 















Figure 4: The branching ratio BR(ii —> bxiff 1 ) as a function of stop (m^) for tan/3=3,5,7,9,ll 
(dashed lines) in mSUGRA model, where ?no=200, m\i 2 = 145, and sign(/r)=+ve. The solid line is 
the leptonic branching ratio for tan/3=3. 



Figure 5: Same as in Fig.4, but for mo=140, m 1 / 2 =180 and sign(/r)=+ve. 
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from collider searches are free from this ambiguity. The upper limit on A 0 can be further 
strengthened from stop searches in the t\ —> b\t channel. 

The conclusion from Fig.5 are similar except for the fact that the allowed values of 
are somewhat larger in this case. From Fig.5 it follows that tan/? ~ 5 the current limits 
from RUN-I are valid. Combining the results of Fig.2 and Fig.5 we find , e.g., |%o| ~ 600 for 
tan/?=7 violates the CDF contraint on the stop mass. 

In order to get a viable signal from 4-body stop decay it is important that the BR 
of decays into electrons and muons be significant. This requirement further restricts the 
interesting region of the parameter space. For a given mi /2 this BR decreases as m o is 
increased. This happens simply because the sleptons get heavier. Of course when m o is 
sufficiently large the BR becomes independent of m o since in this situation only the W - 
exchange diagram contributes. In both Fig.4 and Fig.5, we find that it is possible to obtain 
the leptonic BR > 10% for suitable choices of the SUSY parameters. This may be adequate 
for the detectability of the 4-body stop decays as we shall see in the next section. 


IV. Stop production: Signal and Backgrounds 


In hadron colliders, stop pairs are produced via gluon-gluon fusion and quark-antiquark 
annihilation as, 

. 9 . 9 , qq -> (7) 

The production cross section depends only on the mass of U without any dependence on the 
mixing angle in the stop sector, since it is a pure QCD process. The total pair production 
cross section at the Tevatron for ^/s = 2 TeV is ~ 15-0.3 pb which is 40% larger than the cross 
section for yfs = 1.8 TeV, for the range of rrij A ~ 100-200. The QCD corrections enhance 
this cross section by ~30% over most of SUSY parameter space accessible at Tevatron [|4 . 

We investigate the signal of stop pair production in the channel 7+jets(> 2)+/Z T , assuming 
that one stop decays leptonically and the other hadronically, i.e , 


U -*• bxi^i ; t{ -»• bxiqq', 


( 8 ) 


where q=u,d,c,s and i = e, p. The same event topology also appears from top pair produc¬ 
tion, 


pp —> ££; t —> bqq ', t —> bqq 


( 9 ) 


with a -/y branching ratio suppression. Other main sources of non negligible SM backgrounds 
come from IV boson pair production accompanied by QCD jets, pp —> W + W _ , W+2jets, 
WWj, with one W —> and jets coming either from the other W or from QCD shower. 

We have analyzed the signal and background cross section using parton level Monte Carlo 
simulation without taking into account the fragmentation effects of jets. In our calculation 
we set renormalisation and fragmentation scale, Q 2 = s and used CTEQ4M |25| for the 


parton distributions in proton. The energy of the visible particles in the signal, depends on 
Am = m— rn^o ; the larger the value of Am, the harder is the momentum of leptons and 
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jets in the final state and one obtains a better efficiency of the cuts. We computed the signal 
cross section for a wide range of values varying m^o by fixing r, r = where r> 1 . 

l rn x 3 


The cross section for the process pp —> W+ 2 jets has been estimated using MADGRAPH 
26| . We have cross checked the cross sections with the numbers quoted in Ref. |27|] and 


found that they are consistent within a few percent. Similarly, we also generated the process 
pp —>W W j using the same code. 

For the selection of events, we use the following cuts for the transverse momentum px, 
the rapidity 77 and missing energy j) T . The lepton or jet isolation is selected using a cut on 
A R = \/Acj) 2 + A T] 2 , where Acj) and A// are the difference of azimuthal angle and rapidity 
respectively, between two jets or one jet and one lepton. 


1. Number of jets, n 3 > 2, if they satisfy the pip >20 , \rj 3 \ <2.5 with A R(j,j) > 0.7. We 
require that there be at least one tagged b jet. 

2. Leptons are selected, if pip >10 and \r)g\ <2.5 with A R(£,j) > 0.5. 

3. We require missing energy f> T > 25 . 

4. We require azimuthal angle between lepton and the direction of missing momentum 
be 45° < A <t>{i,i> T ) < 160°. 

5. We demand Ht < 120 , where Ht = pip + i> T - 


Here cuts (1-3) are event selection cuts where as the cuts (4-5) are for background rejec¬ 
tion. We have imposed comparatively harder cuts on jets to minimize the background from 
W boson production. We have noticed that among the isolated jets, almost one is always a 
taggable b jet. 

I 11 Table 1, we show the response of the background processes to these cuts. I 11 background 
processes, the source of lepton and invisible momentum are W decays. Since, in the tt 
background process, particularly, eq|| the final state particles are expected to be more 
energetic relatively to the particles in the signal, the cut [5] kills the tt background by 
~ 30% without much affecting the signal for ~100 . However, it is evident from Table 1, 
the requirement of Ht < 120 is not helpful to suppress the SM backgrounds from W boson 
production with or without jets. The signal cross section, as for example, for mq =120 and 
m^o =60, turns out to be 0.39 pb. Evidently, the level of signal cross section is much below 
the level of background. To reject backgrounds further, we exploited another kinematic 
variable, the transverse mass constructed from the transverse momentum and and missing 
energy [15] : 

m T = sj2pip^ T { 1 — cos[A 0(£,^ t )] (10) 


Fortunately, in top and W backgrounds, the lepton and much of the invisible momentum 
originate from a single W decay leading to a Jacobian peak which is absent in the signal 
process. We studied and found that the mx distribution for the signal is highly populated 
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towards lower values. Therefore, by demanding rrir <45, the background cross section is 
reduced severely. This is reflected in Table. 1, in the the last row. As for example, the most 
dominant background after cuts [1-5] is from W+2 jets. It is suppressed by one order of 
magnitude because of ttit cuts where as tt is suppressed by a factor of ~3 with a modest loss 
in signal cross section. As we see from the last row, the combined background cross section 
turns out to be 0.22 pb leading to ~ 440 events per year with the expected lumonisity 2 fb -1 
per year. 


Cuts 

tt 

W + 2jets 

ww 

WWj 

[1-3] 

.82 

11 

1 

.13 

[4] 

.32 

8 

.72 

.10 

[5] 

.25 

7 

.64 

.09 

m T <45 GeV 

.08 

.06 

.07 

.01 


Table 1: The Background cross sections (in pb) after implementation of kinematical cuts 
described in the text. Respective branching ratio suppression factors are included. 

We have computed the signal cross section for some representative values of r=2, 1.5, 
1.3. In Fig.6, we present the signal cross sections(solid lines) subject to all kinematic cuts 
as described above including the cut on m T . The suppression due to the BR of 4-body 
decay is not included. This will be discussed below in the context of different models. For 
comparison, we also give the stop pair production cross sections(dashed lines) in the same 
figure. Recall that the hardness of the final state particles in the signal depends on Am. 
Larger values of Am mean, higher value of acceptance efficiency. This explains, the rapid 
decrease of signal cross section as shown in Fig.6 for smaller value of r, say 1.3. 

To convert this signal cross section into event rates, we have to multiply the signal cross 
sections in Fig.6 by the respective branching ratio suppression for t\ decay, i.e by 6b r =2. 
Br(ti —> bxi^e)- Br(ti —> bx\qq') =e,/x). In Sec.Ill, we have discussed the decay rates of 
stop in these channels for representative values of SUSY parameters. It follows from Fig.3-5 
that the value of e& r may vary from a few percent to ~ 20%. For the parameter set in Fig.3 
we find that m^o ?»87. Thus for 777^=135, r «1.5. For this mq, e b r — 20%. On the other 
hand from Fig.6 the corresponding signal cross-section is 0.09 pb. Putting all these together 
one can expect for 2 fb -1 of integrated luminosity approximately 36 signal events (S). From 
Table 1 the total number of background events (B) is estimated to be 440, leading to a rather 
modest S /VB ratio. 

To reject the background further, we can impose an additional requirement in event 
selection by tagging one of the jets as a b jet. We know, the major upgradations of CDF and 
DO detectors are underway. Among the improvements the capability to trigger on displaced 
vertices from b quark decays using a precise microvertex detector, is very encouraging. It 
may be possible, to achieve a sufficient b-tagging efficiency, like ~50% or so [|]. Therefore, 
exploiting this facilities we can improve S f\fB by demanding at least one b jet in signal events 
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at the cost of paying a price for b-tagging efficiency in the event number. This requirement 
leads to suppression of the backgrounds from W boson production to an enormous amount 
as there is no b jet in the resulting final state. However, there is a chance that an ordinary 
quark jet may fake as a b jet, but this probability is very very small J28j. The, requirement 


of at least one tagged b jet in the signal will help to get rid of the background process, except 
for the ones from tt, which is not at a negligible level. Nevertheless we gain in S /\fB by 
paying a price for b-tagging efficiency by a factor of 2 in the number of signal events. 

Turning to the mSUGRA parameter set in Fig.4 with m^o ~ 54 we find that for rUOO, 
r ~2. A relatively large signal cross-section as given by Fig.6 leads 100 events in spite of the 
suppression factor eb r — 20%, as can be read off from Fig.4 for tan/3=3. The S /\[B ratio 
thus obtained is quite encouraging even without b tagging. From Fig.5 it is clear that e& r ~ 
20% can be easily realized in the low tan/3 scenario. 

We now discuss the discovery limit of t\ in the channel considered in a fairly model 
independent way. We have seen above that e b r ~ 20% can be realised in a variety of models 
with tan/3=3-4 and this will be used in our analysis as a representative value. It is now 
easy to see from Fig.6 that for 777^=120 one ends up with S /\[B= 4(2) for r ~2(1.5) and an 
integrated luminosity of 2 fb -1 . For r ~2 the search limit can be extended to ~150 for this 


6br • 

The prospect looks even better for the high luminosity option with an integrated lumi¬ 
nosity ~ 15 fb -1 which is expected to accumulate after few years of running. In that option, 
discovery limit of might be better in the proposed channel, even for lower values of r 
and ebr- 

Before ending this section, we note that there is a possibility of detecting the 4-body 
decay of t\ through the jet plus missing energy channel with or without b-tagging. This 
final state occurs in stop pair production when both the stops decay hadronically, eq. ^]. 
This channel looks promising since the corresponding suppression factor eb r is rather mild. 
However in considering this channel, one has to worry about the huge QCD background. 
Nevertheless, it is worthwhile to examine the observability of t\ in this channel. We have not 
discussed this in our present analysis as it is beyond the scope of parton level calculations. 


V. Conclusions 

In conclusion we reiterate that the 4-body decay of t\ BO may indeed play an important 
role in stop searches at the upgraded Tevatron if it happens to be the NLSP. Especially for 
low values of tan/3(=3 - 4) this may be the main discovery channel. 

We have delineated the mSUGRA parameter space where the t\ is the NLSP. For a given 
m (h mi/ 2 , tan/3, sign(/i), a range of values of the trilinear coupling A 0 yields a stop NLSP. 

We have studied the BR of the four body decay of t\ in both mSUGRA and MSSM 
models. We found it to be numerically significant and sometimes even the dominant decay 
mode in a large region of parameter space with low tan/3 (3 ~ tan/3 ~ 5). Here the much 
studied loop induced 2-body decay of t\ is suppressed (Fig.3-5). The 4-body leptonic BR of 
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m~ tl (GeV) 


Figure 6: Cross sections(pb) for stop pair production(dashed line) and signal cross section(solid 
lines) subject to all kinematical cuts, including cut on mx, for some representative values of r, 
where r = —Here the branching ratio suppression due to the stop decays is not included. 


t] which is an essential ingredient of stop search in the channel proposed by us is also found 
to be appreciable. 

In view of the 4-body decay the current limits on rrif A from Tevatron RUN-I based on 
the assumption that the loop induced decay of t\ occurs with 100% BR, needs revision and 
becomes somewhat model dependent. It follows from Fig.4 that the current limit -102 
for o=50 will be significantly relaxed for tan/? ~ 3, since the corresponding BR of the 
loop induced decay is indeed small. For 3- tan/? -5, the 4-body decay is likely to have 
nontrivial impact on stop searches at the upgraded Tevatron. For tan /? ~ 5 the loop decay 
practically overwhelms the 4-body decay. For such values of tan/? the curent limits from 
Tevatron R.UN-I are valid and strong upper limits on % 0 can be placed in mSUGRA models 
for a given set of SUSY parameters as is clearly seen from Fig.l and 2. These limits are more 
restrictive than the ones obtained from the CCB condition or from the requirement that y? 
be the LSP. Thus while the constraint on the mo-mi/ 2 plane from squark-gluino searches ||] 
is fairly independent of Aq, this parameter can be strongly constrained from the negative 
result of stop search for each allowed m 0 -mi /2 pair. 

We studied the viability of discovering the stop in the leptons+jets+^'r channel which 
arise when one of the stops decays leptonically and the other decays hadronically via the 
4-body mode. We have listed the Standard Model backgrounds and the kinematical cuts 
( see the Table ) to suppress them. This discovery limits sensitively depends on 6b r which 
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are model dependent in general and on r. Our studies of the BR in a variety of models in 
Sec. Ill indicate that e& r ~20% is a fairly representative choice. Armed with this information 
we have estimated the discovery limit to be 777^=120(150) for r=1.5(2) for an integrated 
luminosity 2 fb _1 without requiring b-tagging. We have discussed qualitatively how the 
search prospect improves if b-tagging and higher accumulated luminosity after the proposed 
luminosity upgrade are available. The possibility of stop search in the jets+/7r channel 
which has a less severe branching ratio suppression (e^) is also discussed. 

Acknowledgement: 

AD acknowledges the grant of a DST, India Research Project (No. SP/S2/k01/97). SPD’s 
work was supported by a fellowship from Council of Scientific and Industrial Research(CSIR), 
India. MG is grateful to A. Djouadi for drawing attention in this topic and he is also thankful 
to the theory group of CERN where part of this project was carried out. 


References 


[1] For reviews of the Supersymmetry, see e.g. H. E. Haber and G. Kane, Phys. Rep. 117 (1985) 


75; M. Drees and S. Martin, CLTP Report (1995) and hep-ph/9504324 


[2] G. Altarelli, F. Caravaglios, G. Giudice, P. Gambino G. Ridolfi, JHEP 0106 (2001) 018. 

[3] ALEPH Collaboration, Phys. Lett. B469 (1999) 303; DELPHI Collaboration, Phys. Lett. 
B496 (2000) 59; L3 Collaboration, Phys. Lett. B471 (1999) 308; OPAL Collaboration, Phys. 
Lett. B456 (1999) 95; for a summary on stop searches at LEP, see: S. Rosier-Lees et al., 
|licp-pli/9901246| . 


[4] T. Affolder et. al, The CDF collaboration, [iep-ex/0106001 and references therein. 

[5] For a review of mSUGRA and for the physics implications at the taoefTevatron Run II, see: 


S.Abel et al., Report of the “SUGRA working group for 
ph/0003154 and references therein. 


‘RUN-II at the Tevatron”, hep- 


[6] J. Ellis and S. Rudaz, Phys. Lett. B128 (1983) 248; M. Drees and K. Hikasa, Phys. Lett. B252 
(1990) 127. 


[7] DO Collaboration (S. Abachi et ah), Phys. Rev. Lett. 76 (1996) 2222; CDF Collaboration 
(T. Affolder et ah), Phys. Rev. Lett. 84 (2000) 5704; for a summary on stop searches at the 
Tevatron, see: A. Savoy-Navarro for the CDF and DO Collaborations, Report FERMILAB- 
CONF-99-281-E (Nov. 1999). 


[8] K.I. Hikasa and M. Kobayashi, Phys. Rev. D36 (1987) 724. 

[9] J.A. Casas, A. Lleyda and C. Munoz, Nucl. Phys. B471 (1996) 3. 

[10] C. Boehm, A. Djouadi and Y. Mambrini, Phys. Rev. D61 (2000) 095006. 


15 










[ 11 ] 

[ 12 ] 

[13] 

[14] 

[15] 

[16] 
[17] 


[18] 

[19] 


[ 20 ] 

[ 21 ] 


[ 22 ] 

[23] 


[24] 

[25] 

[26] 

[27] 

[28] 


A. Djouadi and Y. Mambrini, Phys. Rev. D63 (2001) 115005. 

W. Porod and T. Wohrmann, Phys. Rev. D55 (1997) 02907. 

A. Datta, M. Guchait and K.K. Jeong, Int. J. Mod. Phys. A14 (1999) 2239. 

W. Porod, Phys. Rev. D59 (1999) 095009. 

H. Baer, J. Sender and X. Tata, Phys. Rev. D50 (1994) 4517. 

R. Demina, J. Lykken, K. Matchev and A. Nomerotski, Phys. Rev. D62 (2000) 035011. 

CDF Collaboration (T. Affolder et al.), Phys. Rev. Lett. 84 (2000) 5273; V. Buescher (for the 
DO collaboration), XXXVI Rencontres de Moriond on Electroweak interactions and Unified 
Theories, 2001; B. Olivier, PhD thesis, Universites Paris VI et VII (April 2001). 

A. Djouadi, M. Guchait and Y. Mambrini, Phys. Rev. D64 (2001) 095014. 

H. Baer, C.H. Chen, M. Drees, F. Paige and X. Tata, Phys. Rev. D59 (1999) 055014 and Phys. 
Rev. Lett. 79 (1997) 986; A. Bartl et ah, Phys. Lett. B435 (1998) 118; A. Djouadi and Y. 
Mambrini, Phys. Lett. B493 (2000) 120; A. Djouadi, Y. Mambrini and M. Muhlleitner, Euro. 
Phys. J. C20 (2001) 563. 

S. P. Das, A. Datta and M.Guchait, in preparation. 

L. Alvarez-Gaume, J. Polchinski and M. Wise, Nucl. Phys. B495 (1983); J.M. Frere, D.R.T. 
Jones and S. Raby, Nucl. Phys. B222 (1983) 11; M. Claudson, L.J. Hall and I. Hinchliffe, Nucl. 
Phys. B228 (1983) 501; C. Kounnas et al , Nucl. Phys. B236 (1984) 438; J.P. Derendinger and 
C.A. Savoy, Nucl. Phys. B237 (1984) 307; M. Drees, M. Gluck and K. Grassie, Phys. Lett. B157 
(1985) 164; J.F. Gunion, H.E. Haber and M. Sher, Nucl. Phys. B306 (1988) 1; H. Komatsu, 
Phys. Lett. B215 (1988) 323. 

Particle Data Group (D. E. Group et. ah) Eur. Phys. J. C15 (2000) 1. 

M. Claudson, L.J. Hall and I. Hinchliffe, Nucl. Phys. B228 (1983) 501; A. Kusenko and P. 
Langacker, Phys. Lett. B391 (1988) 29; A. Kusenko, P. Langacker and G. Segre, Phys. Rev. 
D54 (1996) 5824. 

W. Beenakker, M. Kramer, T. Plehn, M. Spira and P. M. Zerwas, Nucl. Phys. B515 (1998) 3. 


H. L. Lai, J. Houston, S. Kuhlman, F. Olness, J. Owens, D. Soper, W. K. Tung and H. Weerts, 
Phys. Rev. D55 (1997) 1280. 


MADGRAPH by T. Stelzer and W. F. Long, Comp. Phys. co mm n n . 81 (1994) 357; HELAS 
by H. Murayama, I. Watanabe and K. Hagiwara, KEK-91-11 (1992). 


A. Belyaev, E. Boos, L. Dudko and A. Pukhov, 


hep-ph/9511306 


See for example, J. Lykken and K. Matchev, Phys. Rev. D61 (2000) 015001 and references 
therein. 


16 



